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In order to study the growth of wind waves in finite depth we extend Miles' theory to the finite 
depth domain. A depth-dependent wave growth rate is derived from the dispersion relation of the 
wind/water interface. A suitable dimensionless finite depth wave age parameter allows us to plot a 
family of wave growth curves, each family member characterized by the water depth. Two major 
results are that for small wave age, the wave growth rates are comparable to those of deep water 
and for large wave age, a finite-depth wave-age-limited growth is reached, with wave growth rates 
going to zero. The corresponding limiting wave length and limiting phase speed are explicitely 
calculated in the shallow and in the deep water cases. A qualitative agreement with well-known 
empirical results is established and shows the robust consistency of the linear theoretical approach. 
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1. Introduction 

The problem of wind-induced wave growth, whose starting points are the Navier-Stokes equa- 
tions, is a formidable one and several approximations and assumptions were done to fully under- 



stand it. The pioneer works are those of |Jeffreys| ( |1924| |1925| ), Phillips (1957) and |Miles| ( |1957[ ). 
From a mathematical point of view, Miles' theory is based on the dispersion relation of the air-sea 



interface and the related Rayleigh equation ( Rayleigh||1880| |Conte fc M iles 1959) |Drazin &: Reid 
1982). Miles' physical mechanism of wave generation by wind assumes that ocean surface waves 



are generated by a resonance phenomenon. Resonance appears between the wave-induced pressure 
gradient on the inviscid shear airflow and the surface waves. Basically, it occurs when the phase 
velocity of the surface wave equals the flow velocity. The usual model has two spatial dimensions, 
the water is assumed deep, the viscosity is disregarded, the equations of motions are linearized 
and air flow turbulence is largely neglected, except in its role in setting up the logarithmic shear 
flow. The air domain is also considered inviscid and the dynamics equations are linearized around 
a prescribed mean wind velocity. The subsequent Rayleigh equation, depending on a logarithmic 
wind profile, is solved by a combination of analytical and numerical methods. For a good review 



on wind-induced wave growth as well as further progress in quasi linear theory see Miles (1997) 



and Janssen (2004). Alternatively, a review of Jeffreys' theory and some of its applications can be 
found in |Touboul et al.\ ( |2008[ ) . 

Despite these severe simplifications, the Miles' theory brings a sufficiently ideal model to in- 
vestigate the physics of the water wave growth problem. It allows, at least linearly, an analytical 
treatement of the phenomenon. 

Nonetheless, Miles' theory is limited to deep water thus restrictive with regards to winds gen- 
erating nearshore ocean waves or shallow lakes waves. This issue is really challenging the physics 
community and the engineering community. The former because the only way to tackle this basic 
question in interfacial fluid dynamics is via an adequate modelling, and the latter because in coastal 
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engineering works the wave field must be influenced by the water depth. Therefore a theoretical 
extension of Miles's mechanism to finite depth is lacking in the field. 

The situation is very different from the experimental point of view. A lot of investigations were 
done on the subject. Especially the experiments in the shallow Lake George in Australia by Young 



and co-workers in references Young (1997), Young (19991, Young & Verhagen (1996a) and Young 



& Verhagen (1996b). They gave one of the first systematic attempts to understand the physics of 



wave-wind generation in water of finite depth. They provided empirical relationships in terms of 
appropriate dimensionless parameters able to reproduce experimental data. Fetch limited wind- 
waves growth experiments in finite depth water confirmed the existence of asymptotic limits to 
growth of wave energy and spectrum peak frequency. The more important discovery was that at 
short fetches the growth rates tended to be comparable to the deep water limit. At longer fetches, 
the growth rate is lower in shallow water than in deep water, and beyond a limit fetch, the growth 
rate vanishes. 

In this paper we deal with the extension of Miles' theory to the finite depth context. We use 
the Rayleigh equation associated to the air motion with the linear dispersion relation adapting 
the method given by Beji & Nadaoka (2004) for the deep water case. An suitable non-dimensional 
finite depth wave age allow us to derive an explicit expression for the wave growth rate. Evolution 
of the latter against the wave age is plotted. It clearly shows a phenomenon of finite depth limited 
growth. The model is able to derive from theoretical expressions (with a good degree of accuracy) 
important empirical results such as the asymptotic limit where wave growth rate goes to zero. 

The paper is organized as follows. In sections ^ and ( 2.1 1 the linear problem in the water domain 
is introduced after a brief presentation of the model. In section (2.2) the air domain is coupled 
with the water domain and make it possible to solve the linear interface problem. We derive the 
linear dispersion relation of the Miles' theory of wind wave generation in finite depth. In section 
([3]) we introduce dimensionless variables and scalings and we obtain the suitable growth rate. The 
linear dynamic is discussed in section Q with special emphasis on its qualitative agreement with 
empirical relations. Section ^ draws the conlusions and the perspectives. Finally in Appendix |A| 
we give some details of the numerical work behind. 



2. The interface problem. 

We are going to study the stability of an air-water interface. Let the fluid particles be located 
relatively to a fix rectangular Cartesian frame with origin O and axes (x,y,z), where Oz is the 
upward vertical direction. We assume translational symmetry along y and we will only consider 
a sheet of fluid parallel to the xz plane, z — characterizes the interface at rest. The perturbed 
air- water interface will be described by z = rj{x,t). The air occupies the r](x,t) < z < +00 region, 
and the water lies between the bottom located at z = —h and the free surface z = r/(x,t). We 
assume the water as well as the air inviscid and incompressible. The air flow will be described by a 
prescribed mean shear flow, only depending on the vertical coordinate z. We assume the dynamic 
to be linear, and disregard the air turbulence, building a quasi-laminar theory. 

2.1. The water domain 

In the water domaine we consider the Euler equations for finite depth. The horizontal and verti- 
cal velocities of the fluid are u(x, z, t), and w{x, z,t). The continuity equation and the linearized 



equation of motion in the water domain read as (Lighthill (1978)) 



u x +w z = 0, p w u t = -P Xl p w w t = -P g - gp w , (2.1) 

where P(x,z,t) is the pressure, g the gravitational acceleration, p w is the water density and sub- 
scripts in it, w and P denote partial derivatives. The boundary conditions at z = —h and at 
z = r)(x, t) are 

w(-h) = 0, r) t = w(0), (2.2) 
P(x,n,t) = P a (x, V ,t), (2.3) 



where P a is the air pressure evaluated at z — rj. Thus equation (2.3) is the continuity of the pressure 



across the air/water interface. As this is a vital assumption for the growth mechanism, we give it 
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a more pleasant expression. So, let us introduce a reduced pressure defined by 

P(x, z, t) = P(x, z, t) + p w gz - P , 



where Pq is the atmospheric pressure. In terms of ( |2.4| equations ( 2.2 )- ( 2.3 ) read 

0, 



■ w z 



PwU t 



-P x , p w w t 



w(-h) = 0, r) t = w(0), 
~P{x,r],t) = P a (x,rj,t) + p w gr\ - P . 



The linear equations system (2.5)- (2.7) can be solved, assuming normal mode solutions as 
P = V{z) exp (i9), u = U (z) exp (i0) , w — W(z) exp (i9), 77 = rj exp (i9), 



(2.4) 



(2.5) 

(2.6) 
(2.7) 



(2.8) 



with 9 — k(x — ct) where k is the wavenumber, c the phase speed and rjo is a constant. Using 

~ and |2~8| 

w(x, z, t) 



equations ( |2.5| , \2.Q\ , \2.7\ and ( |2.8[ ) we obtain 

-ifccsinh k(z + h) 



sinh kh 



?7oexp (i9), 



kc cosh k(z + h) 
u{x,z,t) = — — 77o exp (i9), 



P(X, 2, t) 



sinh kh 
kp w c 2 cosh k(z + h) 
sinh kh 



770 exp (z0). 



(2.9) 
(2.10) 
(2.11) 



The phase speed c is unknown in equations ( 2.9 )-( 2.11 ). To determine c we have to consider the 

(2.12) 



boundary conditions (2.7) (not yet used) and (2.6) which yields 

c 2 kp w r]o exp (i9) coth kh — p w gt]o exp (i9) + Po = P a (x, 77, t). 
In the single-domain problem P a (x, r),t) — Pq and ( 2.12| ) gives the usual expression for c, 

c 2 = cl = y tanh(M). (2.13) 
It is not the case here, and we must use the air pressure evaluated at z = 77 to determine c. 

2.2. TTie air domain 

Let us consider the linearized (inviscid) governing equation of a steady air flow, with a prescribed 
mean horizontal velocity Viz) depending on the vertical coordinate z. We are going to study 
perturbations to the mean flow U(z): u a (x, z,t), w a (x,z,t) and P a {x 1 z,t) (subscript a stands for 
air). So with ~P a (x,z,t) — P a {x,z,t) + p a gz — P Q , p a the air density, and U' — dU(z)/dz we have 
the following equations 

U a ,x + w a , z = 0, (2.14) 
Pa[u a ,t + U(z)u atX + U' \z)w a ] = -P , B , (2.15) 
Pa[w a ,t + U(z)w a ,x] = -P<j, z , ( 2 -16) 

which must be completed with some appropriate boundary conditions. The first one is the kinematic 
boundary condition for air, evaluated at the aerodynamic sea surface roughness zq placed just above 
the interface. It reads 

rit + U{z )r] x = w a (z ). (2-17) 

We choose U(z) to be the logarithmic wind profile. The log wind profile is commonly used to 
describe the vertical distribution of the horizontal mean wind speed within the lowest portion of 
the air-side of the marine boundary layer ( Garratt et aL|[l996 l 



U(z) = U 1 ln(z/z ), Ui = —, KRiO.41, (2.18) 

K 

where is the friction velocitjjf] and k the Von Karman constant. So, eq. (2.17) can be reduced to 

Vt = w a (z ). (2.19) 

f Scaling the wind speed with the friction velocity is not a keystone of this work, and the reference speed 
can be changed. 
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This equation describes the influence of the surface perturbation on the vertical perturbed wind 
speed. Next we assume P = V a {z) exp (iO), u a — 1A a {z) exp (i6), w a — W a {z) exp (i6) and we add 
the following boundary conditions on W a and V a , 



lim {W a + kW a ) = 0, 

z— >-+oo 



(2.20) 
(2.21) 
(2.22) 

that is, the disturbance vanishes at infinity, and the vertical component of the wind speed is 
enforced by the wave movement at the sea surface. Then, using equations ( 2.14 1-( 2.16 1 and (2.221 
we obtain 



lim W a = W , 

z->z 

lim T a = 0, 



w a {x,z,t) = W o exp(i0), 
i 

U a (x,Z,t) = -W a .zGXp(ii 
K 



/>oo 

P a (x,z,t) = ikp a exp(i9) / [U(z') - c]W a (z')dz'. 

J Z 



(2.23) 
(2.24) 

(2.25) 



Removing the pressure from the Euler equations, we find the well-known Rayleigh equation (inviscid 
Orr-Sommcrfcld) 



(U - c)(W a ' - k 2 W a ) - U"W a =0 Vz \ z < z < +oo, 



(2.26) 



which is singular in z c — ZQe CK ^ w " > z$ > 0, where U(z c ) — c. In equations ( 2.23 )-( 2.26) neither 
yV a (z) nor c are known. In order to find c, we have to calculate P a (x,i],t). We obtain 

/>oo 

P a {x,T ll t)=P -p a gr ] + ikp a exp(i9) [U (z) - c]W a {z)dz , (2.27) 

J z 

where the lower integration bound is taken at the roughness height zq instead of z = r\ since we 
are studying the linear problem. Finally, using equation (2.19) to eliminate the term ikp a exp(i8) 

-sk 2 



the equation (2.27) in (2.12) yields 
sk 2 



5(1 



U(z)W a (z)dz - c 2 { 



Wn 



W a {z)dz + fccoth(fe/i)} = 0, 



2.28) 



where s — p a /Pw This is the dispersion relation of the problem. If h — > oo we obtain the expression 
(3.7) found by Beji & Nadaoka (2004 1. The parameter s is small (p a /pw ~ 10 -3 ) and (2.28) may 
be approximated as 

c = c + s Cl +O{s 2 ). (2.29) 

The explicit expression of c\ is calculated in the next section. Therefore, we can find W Q (z) by 
solving (2.26) with c substituted by cq, that is to say, of order zero in s. The method is shortly 
described in appendix [A] 



3. Wave growth rate 

The function W a (z) is complex and consequently c also. Its imaginary part gives the growth rate 
of rj(x, t) defined by 

7 = fclm(c), (3.1) 
where Im(c) is the imaginary part of c. The theoretical and numerical results concerning the growth 



rate -, are studied and computed in terms of two dimensionless parameters S (see Young fc Verhagen 
(1996a) and Young & Verhagen (1996b)) and 9d w defined by 



x 9 h 



(3.2) 



The dimensionless parameter 5 is the more important novelty in relation to the set of non- 
dimensional parameters governing 7 in deep water. This parameter measures the influence of the 



5 



finite fluid depth on the rate of growth of rj(x, t). The parameter 9d w is nothing more than the deep 
water wave age. Now a finite depth wave age 9f d can be introduced as 



1 19 



->/,! = jj-y -\/t&nh(kh) = OdwT 1 



/2 



(3.3) 



where T = tanh(gj-). The form (3.3) defines a depth weighted wave age such that: for a finite 

and constant 9d w we have dr d ~ @dw if S — > oo and Qj d ~ 5 1 / 2 = y/gh/U\ if 5 — >■ 0. In order to 
obtain the growth rate, we introduce the following non-dimensional variables and scalings (hats 
mean dimensionless quantities) 



U = U X U, W a = W W a , Z 



— , c = U\c, t = — t. 
k g 



(3.4) 



Using (3.2) and (3.4) in equation (2.28) and retaining only the terms of order one in s we obtain c, 

s 



'd a: ) 



e dw T^ - °-e dw T x / 2 + °-{Th - e dw T^i 2 }, 



(3.5) 



and with e 7 * = e k " xi ^ t = e ( c )*/ 9 < i «», we have the dimensionless growth rate 7 = ^7 as (dropping 
the hats) 



_ s Tlm(h) T 3 / 2 Im(/ 2 ) 

~ V dw V d w 



where the integrals are defined as follow 



h = 



UW a &Z, h = 



Wadz. 



(3.6) 



(3.7) 



So, we can compute it for a given (5, 9 dw ) set. The 5 parameter does not appear explicitly allowing 
us to indeed compute 7 for an infinite depth, where we have just T — > 1. This gives back Miles' 
theory. 



4. Discussion 



The existence of a finite depth h transforms the unique curve of wave growth rate in deep water 
in a family of curves indexed by S — gh/Ui, i.e., a curve for each value of S. In Figure [l] is shown 
a family of five values of S against the finite-depth wave age. The limit 5 — > 00, corresponding to 
the envelope of the family of curves, is included as well. 

Figure [T] clearly shows that at small wave age the growth rate 7 is almost equal for all values of 
S, the limit being the deep water case. 

As the wave age increases, the finite-depth effects begin to appear. The growth rate becomes 
lower than in the deep water limit, for each value of 6. The growth rates are scaled with 5: for 
a given wave age, the bigger the 8 the larger the 7. Each <5-curve approaches its own (idealized) 
finite-depth wave-age-limited growth as 7 goes to zero. At this stage the wave reaches a final state 
of linear progressive wave with zero growth. In others words, for a given 5 the surface wave does 
not grow old anymore beyond a determined wave age. 

We note that the theoretical curves in Figure [T] are, mutatis mutandi, in good agreement with 
the empirical curves of the fractional wave energy increase per radian T as a function of the inverse 
wave age Uio/Cp in Young (1997) ( see Figs. 3a, 36 3c and 3d). 

In Figure [2| are plotted, against 8, the critical values of the wave age 8j d for which the growth 
rate 7 goes to zero. They obey the relation 



7 fd 



= 5 1 ' 2 . 



(4.1) 



The above relation, found numerically, is coherent with the parameter formulation ( |3.3[ ). It is 
indeed a limiting value for the wave age, uniquely determined by the water depth. This result is 



in excellent agreement with the one given by Young (1997) for Y (from Donelan et al. (1992)) The 



authors have shown from an empirical relationship (formule (6) in reference above) that T goes to 
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zero for 



Cp 



= 0. 



gh_s 

2 / 



0.45 



■u 



(4.2) 



10 



In Figure [2] are also represented data from Donelan et al. ( 2006 ) , from the Australian Shallow 



Water Experiment, wich will be referred to as AUSWEX. A basic fit is also plotted to show the 
general trend. The raw data consists in the water depth h in metres, the friction velocity u*, the 
10 metres wind velocity Uio and the ratio of the former with the measured phase speed c p , U\q/c p . 
For example, it* = 0.44 m.s -1 , h — 0.32 m, Uxo = 11.9 m.s -1 and Uxo/c p — 7.2 gives 5 = 2.7 and 



Ofd — 1.55, which gives a small relative error regarding (4.1 ). All the points give (S,9fd) coordinates 
really close to the theoretical limit. 

With c j d we can calculate the corresponding critical wave length A c . Using (4.1| in (3.3 1 we 
obtain 



= tanh (^-) 



(4.3) 



Relation (4.3| means the wave has entered the shallow water region. In such a limit the range 
of 5/6 dw is: < S/0 2 dw < f ( |Fenton| ( |1979| , |Francius k Kharif] ( |2006[ )). As a result we obtain 
A > A c = 8h. For these values of A, the phase velocity is in the long wave limit i.e., c = \fgh. 
Consequently, if A > A c the wave feels the bottom, the amplitude does not grow anymore, the 
resonance wind/phase speed ceases, and the wave reaches its utmost state as a progressive plane 
wave. 

An analogous phenomenon appears in the deep water limit <5 — > oo where tanh (gr-) is in [0.99; 1] 
for 6/9 dw > 7r, so we have 



_ s . Im(Jx) 
7 _ 2 *■ 2 

z °dw 



Im(/ 2 ) 



}• 



(4.4) 



The A c is in this case A c = 2h. For A < A c , 7 goes to and the sea reaches the final state for which 
we can obtain from (4.4) the phase velocity 9d w 



J dw 




Im(Ji) MJ™ U(z)W a (z)dz) 



Im(J 2 



(4.5) 



In contrast to the usual analysis in wind-induced wave growth (in deep or finite depth), our 
results concern the dimensionless growth rate 7 in terms of the finite depth wave age Ofd instead 
of a finite depth fetch Xfd- Therefore a transformation rule is needed. It can be easily derived from 
the empirical relation introduced by Hasselmann ( Hasselmann et oT]|1973 ) 

D 



f P U- 



10 



9 



3.5 



xg_ 

U 2 
u io 



D 



-0.33, 



(4.6) 



with f p the peak frequency, Uio t ne wind speed at 10 m and x the fetch in metres. Now using 
f p = g/{2TTc Q ), 5 = and the formula in (2.18) we obtain a transformation between 8 fd and the 
dimensionless fetch Xfd 



1 



In 



3 z io 



1+2D 



acU 2 K 2 

10 m and the Charnock constant ac 



(Xfd 



Xfd 



xg^ 

ur 



(4.7) 



0.018 used to determine the roughness length 



with z%o 

z o = ®cu 2 /g. This relation was proposed by Charnock ( 1955 ), and is used since, although other pa- 
rameterizations of the roughness length exist. Johnson's relation, for example, states a dependence 
between Zq and the wave age ( Johnson et aL|19"98 ). Identically we have the following transformation 
rule between the Miles' dimensionless growth rate f3 and dimensionless 7 



27, 
s 



,3 rpl/2 



(4.8) 



where we took f3 as it is usually defined, with the dimensions, Im(c) = Co§/3(— ) . This is a 
straightforward definition of Miles' f3 in finite depth. Its evolution is shown clearly in Figure [3j 
showing the usual deep water trends, and the new finite depth limits. The effects of depth are 



7 



Growth rate vs. wave age for several 5s and deep water 




2 4 6 8 10 
Wave age 0fd 



Figure 1. Evolution of the growth rate in semi-logarithmic scale. The curves that goes to zero almost 
straightly are the finite depth ones. From left to right, they match 8 — 4,9, 25, 49, 81. We see that for each 
depth, there is a wave-age limited growth. The true deep water limit, also computed, is approached for 
small Ofd and matches Miles' results. 

critical. The /3 is almost constant for small 9fd, as usual, but it goes to zero extremely fast when 
the finite-depth wave-age limit is close. 



5. Conclusions and Perspectives 

In this work we built a theoretical extension to the finite depth domain of Miles' theory of 



wave generation by wind (Miles (1957)). In order to compute wave growth rates in finite depth 



we adapted the method introduced by Beji & Nadaoka (2004) for deep water. We defined a finite 
depth wave age useful to study families of wave growth rates, each member of the family tagged 
by a water depth. 

From a numerical study based on purely theoretical relationships it is shown that wave growth 
rates go to zero for long waves in shallow water. For small wave age the wave growth rates behave 
as in deep water, regardless of the actual depth. These results are the qualitative analogous of 
the well known empirical (or semi-empirical) results concerning wave energy increase per radian in 
function of the inverse wave age (Young (1997) Young ( 1999| ) ). We have shown, for the first time, 
curves of the original Miles' dimensionless growth rate j3 in finite depth. This growth rate drops 
from a deep water behaviour to zero extremely quickly as the wave age reaches its limit. 

Miles' theory studies exponential growth of the wave amplitude of a linear and uniform monochro- 
matic wave train. Nowadays nonlinear modulational analysis beyond Miles' monochromatic theory 
were carried out. In reference Kharif et al. (2010), the authors report the behaviour of the Ben- 
jamin Feir instability when dissipation and wind input are both taken into account. Within the 
framework of weakly nonlinear modulated wave train, with this same theory as the foundations, 
we are going to study the influence of the wind in the finite depth Benjamin- Feir instability (for 
kh > 1.363). Inclusion of small bottom friction and constant vorticity are also under consideration. 
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Curve of low- valued 7s 



7 fd 




Computed limit o 

Theoretical limit 

AUSWEX data o 
AUSWEX fit 



FIGURE 2. Parameter curves where the growth rate is alm ost zero. The theoretical limit is (4.1 1. The 
AUSWEX data was taken from |Donelan et al\ |200oj), and the AUSWEX fit is 6 fd = OM^ 2 . 



Miles' (3 vs. wave age for several 5s and deep water 
4 




8 = 4 
8 = 9 
8 = 25 
8 = 49 
8 = 81 
S — > 00 
Miles 



Figure 3. Evolution of Miles' (3. Every curve is calculated with the same z.o, the differences are only due 
to the various depths. As we see, the finite-depth effect is critical, and for high 8, we reach a deep water 
behaviour. 



Appendix A. Rayleigh equation 

We recall the Rayleigh equation 

(U - c)(W^ - k 2 W a ) - U"W a = Vz \ z a < z < +00 



(Al) 



which is singular in z c > z > 0, where U(z c ) = c. This equation underlies the most essential 
mechanism of flow stability. This equation is singular only for the zeroth order in s, where c = 
Co + o(s), which is real. The value of ci, hopefully complex, is then found with the dispersion 
relation (2.28). Nevertheless, this is computational, and do not change the fact that we search 



a complex c eigenvalue in (Al). We have to prescribe a flow U(z) allowing instability, or more 
precisely, not forbidding it. Rayleigh's inflexion point theorem, and subsequent Fjortoft theorem 



states that U{z) is bound to have one inflexion point to, at least, not forbid instability (see Fjortoft 



( 1950 1 ) . Usually, for this theorem, the domain of W a is [0, +00 [, and the boundary conditions are 



the vanishing at each bound. But our lower bound is zq, which is nonzero, and we have a forcing 
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of W a in zq. Taking this into account, we derive the following constraint 

'■+ 00 U"{z) 



Im(c) 



\U{z) 



\W a (z)\ 2 dz = -Im( lim W*(z)W a (z)). 



(A 2) 



Where we see that if zq — > and W a vanishes smoothly at the boundaries, considering that U"(z) 
is monotonous, indeed the r.h.s vanishes and Im(c) must be zero. Here, it can be nonzero. So, the 
condition found allows an exponential growth of the free surface 77(2;, t), like a mechanical oscillating 
system forced into one of its normal modes ( Conte fc Miles||1959[ ). We use now a semi-numerical 
recipe to solve ( A 1 ) for W a following the method introduced by Beji and Nadaoka. We first develop 



(All in the z c -neighbourhood, assuming 

U{z)^U'(z c )(z-z c ) + c, U"{z)nU"(z c ), (fc 



U'(z c ) a ^ Q 



(A3) 



U"(z c )' 

These are fairly true in the logarithmic wind profile case. After some algebra, these assumptions 



transform (A 1 1 in a Bessel equation of order one, whose solutions are known to be a linear combina- 



tion of the two first-order Bessel functions. So, with the weighted-centered, dimensionless variable 
z n = - Viz^^H = - - 1 we find 



U'( Ze ) 



W a (z p ) = v ^(C 1 J 1 (2y^) + C a yi(2V^)) . 



(A 4) 



Where C\ and C 2 are complex constants. Here, the function W a is z p -dependent. The numerical 
solution that we seek can be written as 



W a = C l W l +C 2 W 2l 



(A 5) 



where W\ and W 2 are unknown and independent. W\ (resp. W 2 ) is found by integrating ( A 1 ) with 
Ji (resp. Yx) for initial function value and slope value at z c . The first thing we need is to avoid 
the numerical singularity. We achieve this by introducing a small parameter e. We choose the two 
initial points to be z^r — (1 ± e)z c , and evaluate the Bessel functions in these points. As we notice, 



1/2 ! 

the function z p Y\{2zp^) becomes complex, with a negative imaginary part, when z < z c . Its 
derivative becomes complex, with a positive imaginary part. 

We have to take the complex-conjugates of the initial value and initial slope value of these 
functions to get a positive growth ratc[f] (Drazin & Rcid 1982). Therefore, after integration, we 
get Wi and W 2 . But this doesn't take into account integration constants. These are set using the 



1/2 



boundary conditions of surface forcing and vanishing at infinity. So, using ( A 5 1 in ( 2.20 ) and ( 2.21 ) , 



we obtain a simple algebraic system which allows us to determine C\ and C 2 . The "infinite" value, 
for computation, is the one from which the constants are stable enough. The relative error on the 
constants is proportional to the error on the growth rate. Here, a relative error of 10 -5 is taken. 



As a result, we obtain W a (z), and we can evaluate the integrals in (3.6) for any parameters set. 
The results are shown in Figure [T] 



f Taking the non-conjugate values will simply give a minus sign to Im(c). 
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